Non-locality of Foldy-Wouthuysen and related transformations for the Dirac equation 
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Non-localities of Foldy-Wouthuysen and related transformations, which are used to separate posi- 
tive and negative energy states in the Dirac equation, are investigated. Second moments of functional 
kernels generated by the transformations are calculated, the transformed functions and their vari- 
ances are computed. It is shown that all the transformed quantities are smeared in the coordinate 
space by the amount comparable to the Compton wavelength A c = h/mc. 
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I. INTRODUCTION 

The Dirac equation for relativistic electrons, in spite of 
its fundamental importance in physics, is far from being 
thoroughly understood [l[. One of its special features is 
that, even in the absence of fields, there exists a spec- 
trum of positive and negative electron energies. This 
peculiarity is a source of various problems. One of them 
is a phenomenon of Zitterbewegung (trembling motion) 
which arises from an interference of positive and negative 
energy states There exist attempts to circumvent 

this duality. One of them is the Foldy and Wouthuysen 
transformation (FWT) which, for the case of no exter- 
nal fields, allows one to break the Dirac equation into 
separate equations for positive and negative energies @. 
In their original paper, FW remarked that a functional 
kernel, which transforms functions from the original rep- 
resentation to the FW representation, is characterized 
by a non-locality in coordinate space of the order of the 
Compton wavelength A c = h/mc. Rose in his book Q 
put this statement on a quantitative basis by showing 
that the second moment of the kernel is equal to (3/4)Aj! 
(see below). This result was not followed by other in- 
vestigations and it is by now not well known. Also, the 
moments are often used in statistical physics but they 
are a rather unorthodox way to characterize quantum 
mechanical properties. It was shown later that the FW 
transformation is not unique. In other words, there ex- 
ist other transformations capable of separating positive 
and negative energy states in the Dirac equation, both in 
the absence of fields and in the presence of an external 
magnetic field [§-[llj]. I n the following we consider, in ad- 
dition to the FWT, a two-step transformation proposed 
by Moss and Okninski (MO, Ref. pf). 

Let us consider, as a matter of example, an average 
value of the time-dependent velocity operator v(t) in the 
Heisenberg picture, 



($|fi(t)|$) = ($U f \Uv(t)U^\U$), 



(1) 
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where v(t) = e l H D t/h- e -iH D t/h ^ jj D ig the Dj rac Hamil- 
tonian and |<f>) is an arbitrary Dirac spinor. Accord- 
ing to the Dirac equation one has Vj = cay, where &j 
(j = x,y,z) are the standard Dirac 4x4 matrices. As 
is well known, aj do not commute with Hp, so the 
velocity depends on time also in the absence of fields, 
which results in the Zitterbewegung mentioned above. 
The operator U in Eq. (fTJ) stands for a unitary transfor- 
mation, for example the FWT. Equation (JlJ expresses 
the well known fact that a unitary transformation does 
not change an average value. In our case it means that 
a transformation will not eliminate the Zitterbewegung, 
which makes physical sense. We wrote down Eq. (JTJ to 
illustrate symbolically that, if one transforms the oper- 
ators: v(t) = Uv(t)W, one should also transform the 
functions: $ = U$>. Important operators transformed 
according to FWT are given in the original FW paper 
and quoted in many textbooks. Here, we concentrate on 
the transformed functions which were not analyzed in the 
literature. 

Our work has two objectives. The first is to char- 
acterize in various ways the non-locality of the Foldy- 
Wouthuysen transformation for the Dirac equation. In 
doing this we continue the work of Rose but use means 
more typical for quantum mechanics. The second objec- 
tive is to investigate, as a matter of example, another 
"separating" transformation and to compare its proper- 
ties to those of FWT. Such a comparison should give 
an idea of what one can expect of various " separating " 
transformations. 

Our subject is of relevance for two reasons. First, 
with "the rise of graphene" there is nowadays a great 
deal of interest in relativistic-type wave equations fl3l — 
[l5j . Second, it is now possible to simulate the Dirac 
equation with the trapped ions and cold atoms interact- 
ing with the laser radiation, where one can tailor much 
more "user friendly" values of the basic parameters mc 2 
and c 

pun. 

In fact, a proof-of-principle experiment 
simulating 1+1 Dirac equation and the result ing Zitter- 
bewegung was carried out by Gerritsma et al. [201 ] . If one 
deals with effective parameters by employing narrow-gap 
semiconductors, the non- locality is determined by an "ef- 
fective" Compton wavelength Xz = h/m*u, where m* is 



2 



the effective electron 
velocity. In narrow-g 
m* ~ 10~ 2 m and u 
can be a few orders 
vacuum, see Ref. [if 
corresponds to the ' 
a similar maximum 
be even larger. 



mass and u is the maximum electron 
;ap semiconductors one can deal with 
~ 10 8 cm/s, so that the non-locality 
of magnitude larger than that in a 
j|. The band structure of graphene 
ultra-relativistic" case (no gap) and 
velocity u, so the non- locality may 



II. TRANSFORMATION KERNELS AND 
NON-LOCALITY 

In this section we establish forms of Foldy-Wouthuysen 
and Moss-Okninski transformations and define corre- 
sponding transformation kernels. Next we characterize 
nonlocal properties of the kernels by calculating their mo- 
ments. The Dirac equation for a free relativistic electron 
is Hp = cap + (3mc 2 , where p is the momentum and ay 
and P are 4x4 Dirac matrices in the standard notation. 
The eigenenergies are ±E p with E p = \J m 2 c 4 + c 2 p 2 . 

For the field-free case the momentum components are 
good quantum numbers p — > p. In the FW transfor- 
mation an initial function |$) is transformed to the new 
representation |$) with the use of the operator Ufw de- 
fined as 



Ufw = 



E p + (3H P 



^/2E p (E p + mc 2 ) 



(2) 



There is |$) = Ufw\$)- To find the MO transformation, 
we follow Moss and Okninski .12] and introduce first the 
unitary operator 



(3) 



which transforms the Dirac equation to the form H' p = 
cap + Smc 2 , in which 



S = a x d y a z (3 



-i 

-i 

1 
i 



(4) 



The Hamiltonian H' p has zeroes on the diagonal. For 
a four-component function \9) there is |W) = V^). 
The two-step MO transformation is defined using the H' p 
Hamiltonian in the following way 



(5) 



In the field-free case the operators U fw and Umo are 
number matrices. Following Foldy and Wouthuysen we 
introduce a function transformation from the old to the 
new representation (r\<t) = (r\UQ). Inserting the unity 



operators we obtain 

d 3 pd 3 p'd 3 r' 



$(r) 



(27rft) 6 
d 3 pd 3 r' 



r\p)(p\U\p')(p'\r')(r'\<£) 



(27rft) 3 

which defines the transformation kernel 

ip(r-r')/h 



K(r,r') 



(2tt^) 3 



-U p d 3 p. 



(6) 



(7) 



In Eqs. ([6]) and (0 U stands for each of the two operators 
defined above and U p = (p\U\p')S(p — p'). The kernels 
are 4x4 matrices. 

As mentioned above, the transformation kernels have 
a nonlocal character. We also remarked that Rose Q 
described this feature quantitatively by calculating the 
second moment of the FW kernel and showing that its 
smearing is given by the Compton wavelength A c . Below 
we evoke his calculations, as it is not easily accessible 
by now, and compare his result with the corresponding 
quantities for the MO transformation. The moments are 
rarely used to describe quantum mechanical behavior, for 
this reason we do not limit our subsequent considerations 
to the moments but also consider other quantities. 

All the elements of matrices defined in Eq. ([7]) are 
nonlocal functions of R = r — r'. To estimate the 
non-locality we calculate the zeroth and second mo- 
ments MW of Kpw(R) and K M o{R) matrices 



= I Kij(R)R n d 3 R. 



(8) 



Let us consider first the (1, 1) element of Kfw(R)- For 
the zeroth moment there is 



M 



(0) 

li 



d 6 R 



^ipR/h 



E p + mc 



(27rft) 3 ^2E p (E p + mc 2 ) 



d 3 p. (9) 



Changing the order of integration one has J e tpR ' ti d 3 R — 
{2ttH) 3 8(p). Then the integration over d 3 p is trivial and 

one obtains = 1. The second moment of the (1,1) 

element is 



Mi? 



d R 



jpR/h 



R 2 d 3 p. 



(27rft)3 y/2E p (E p 
Integrating twice by parts one obtains 

£ ipR/h 



M? 



d 3 R 



(27rft) 3 
E n A 



h 2 (E p 



y/2E p {E p + mc 2 ) 



(10) 



d 3 p 



V / 2E P {E P + mc 2 ) 



p=0 



l^l- (11) 
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In the same way one calculates moments of the remaining where 



matrix elements. We obtain M 



(0) 
FW 



I and M 



(2) 
FW 



(3/4)A 2 /, where I is the identity matrix. These are the 
results of Rose Q ■ 

Thus the Kfw kernel is smeared in the R = r — r' 
space by the amount of A c /2 in each direction. The ker- 
nel can be considered to be a transformed Dirac delta 
function, cf. Eq. ([5]). Hence, one can subscribe to the 
statement of Foldy and Wouthuysen that "a wave func- 
tion which in the old representation corresponds to a 
state in which the particle was localized at one point, 
passes over in the new representation into a wave func- 
tion which corresponds to the particle being spread out 
over a finite region" . The physical meaning of this result 
is discussed in Refs. [1, 0|- 

We carry similar calculations for the Moss-Okninski 
transformation and calculate the zeroth moment of 
the (1,3) element of Umo 



M 



(0) 



d 3 R 



CPz 



(2^)3 y/2E p 



d 3 p- 



i 

71' 



(12) 



For the second moment of (1,3) element of Umo we in- 
tegrate twice by parts and obtain 



M 



(2) 
13 



<TR 



= -ti 2 V 2 



e ipR/H 

(27rfr) 3 



CPz 



V2E, 



d 3 p 



p=0 



p 

3*A e 2 
V2 



(13) 



Calculating the remaining moments one finally obtains 

S), (14) 

S. 



M MO 



A 2 



V2 



(15) 



Thus the kernels of both FW and MO transformations 
have finite non- localities of the order of A c , but they differ 
somewhat from each other. In this sense the FW trans- 
formation is somewhat more "compact" than the MO 
transformation. 



III. TRANSFORMATION OF FUNCTIONS 

We further investigate properties of FW and MO trans- 
formations by studying the transformed functions. Let 
us consider an initial wave function ^ — /(r)(l, 0, 0, 0) T , 
where f(r) is normalized. Then the transformed function 
for the MOT is #'(r) = %(r) + %(r), in which 



1 



(16) 



1 f d 3 p e ipr l h f p 
2j (2^)3 E p 



mc + icp z 

"V ' I ( | 7 , 

imc 2 + cp z ' ' ' 



cp 



+ 



- lpr '/ h f{r')d 3 r'. 



(18) 



The prime in ^'(r) means that we transform the function 
#'(r) = VV(r), see Eq. Q. 

If the initial function $(r) = /(r)(l, 0, 0, 0) T is trans- 
formed according to FWT, there is 



$(r) = 



d 3 p e l P r / h f p 



E p + mc 2 


-CPz 

—cp + 



(19) 



where N p = y / 2E p (E p + mc 2 ). The asymmetry between 
two upper and two lower components of $(r) arises from 
the asymmetry of components in the initial wave packet. 

Let us consider first the initial function in form of the 
delta function: *(r) = 5(r)(l, 0, 0, 0) T . Then f p = 1 
and the transformed function is obtained from Eqs. (|16p 
and (|17p in terms of four integrals Do , D x , D y , D z 



The first integral is 




(20) 



D 



d 3 p 



Jpr/h 



2{2nK) 3 E p 



(2tt) 2 A3 r 



Ki(r/\ c ), 



(21) 

where Kx(z) is the modified Bessel (MacDonald) func- 
tion. To get Eq. $ZT$ we used identities (|AT|) and (|A2|) 
in Appendix. The second integral is 



D, 



d 3 p 



2{2-Khf 



E n 



(22) 



This integral is divergent, so we separate it into the di- 
vergent and convergent parts and obtain 



D z = 



d?pe ipr ' n , . 

w sgn(ft) 

d 3 pe i P r / h cp, 
2(2irh) 3 



\J m 2 c 4 + c 2 p 2 



sgn(p z ) 



2X r 



l -8{p)+A{r/\ c ), 



(23) 



where p — (x, y). We marked the integral in the second 
line of Eq. (|23|) by A(r/X c ). This integral is convergent 
since the integrand has no singularities and for large p it 
decreases as p~ 2 . The integral A(r/X c ) can be expressed 
in terms of the Anger functions and it decays exponen- 
tially with \r\ with a characteristic length A c . The inte- 
grals D x and D y can be obtained from D z replacing z 
by x and y. In conclusion, it can be seen that the delta 



4 



function subjected to the MO transformation becomes a 
function having a finite width of the order of A c . 

Similar calculations can be performed for a function 
$(r) = £(r)(l, 0, 0, 0) T subjected to the FW transforma- 
tion. The transformed function is a combination of four 
integrals: I = (B o ,0,-B z ,-B x - iB y ) T , see Eq. (|T9| . 
Setting f p =l and separating out the divergent part we 
obtain 



So 



d 3 pe tpr/h 

6(r) f d 3 k 



V E p {E p + 



1 



1 



f G(k), 



V2 J V2(2tt) 3 yiTAffc 
where G(k) = ( 1 + Jl + + X 2 c k 2 



(24) 



We 



changed p to Hk. The function G(k) varies slowly from 
G(0) = 1/(V2 + 1) 0.414 to G(oo) = 1/2. Thus we 
may approximate _Bo by setting G{k) to have a constant 
value Cg. This gives 



B 



Co 



A 



J(r) 

V2 V2(2^) 2 A3 r 



£ ^i(r/A c ). 



(25) 



This final result is very similar to that for the MOT, see 
Eqs. ([201) and (JUJ). The integrals B x , B y and B z can 
be calculated in a similar way to that described for the 
MO transformation. Thus for both MOT and FWT the 
initial delta function transforms into functions of finite 
width of the order of A c . 

In the above analysis with the initial delta function 
the interpretation of results is somewhat difficult because 
of divergent and singular final integrals. To avoid these 
problems we take the initial function in the form of a 
Gaussian packet having a finite width d 

exp(-r 2 /2d 2 ) 



f(r) 



^3/4^3/2 



(26) 



The packet is normalized according to J \f(r)\ 2 d 3 r = 1. 
Now there is f p = (2dy/n) 3 / 2 exp(-p 2 cZ 2 /2?i 2 ) and all the 
integrals in Eqs. (| 1 T|) and (fl~9|) are quickly convergent. 
The transformed function is again expressed in terms of 
four integrals To, T x , T y and T z 



f(r) 









- iT z 






iT x 


-Ty 


-* 


+ 


iT 


+ T Z 


o / 


\T X - 


-iTy 



(27) 



The first integral is [see Eq. ([IT)] 



To = V2mc 2 {dV^f /2 J 



d 3 p 



(2Trh) 3 ^ m ^c 4 + c 2 p 2 



(28) 



Integrating over the angular variables in the spherical 
coordinates we obtain 



To = 



(cV^F) 3 / 2 f°° sin(fcr) e^^ ' 2 k 2 dk 



kr 



yr+FA 2 



(29) 



Applying the identity: 1/a = J x oo exp(~ar] 2 )dri/^/Tr, in- 
tegrating over k with the use of formula (|A3[) . we get 



To = 



3/2 



V2vr 2 



e -V 2 e -r 2 /(2d 2 +^ 2 X 2 c ) 

4(d 2 /2 + ?7 2 A 2 ) 3 / 2 



drj. (30) 



The presence of e~ v term in Eq. (|30l) reduces the range 
of integration to |?7| < 3. In the limit of d 3> A c we can 
neglect ?7 2 A 2 as compared to d 2 and obtain 



To 



(dV5F) 



3/2 -r 2 /2d 2 



V2ir 2 4(d 2 /2) 3 / 2 7- 



e- r > dr) 



e -r 2 /2d 2 
27T 3 / 4 d 3 / 2 ' 



(31) 

Thus, there is no widening of the wave packet in the limit 
of large widths d, and the transformed packet is almost 
identical to the initial one. In the opposite limit of very 
narrow packets: d <C A c , we may neglect in Eq. (|3T)]) the d 2 
term under the integral sign. After the substitution \ = 
rj 2 one gets 

r^Wggr^^x. (32) 



V2vr 2 A 3 



A X 2 



Using identities (|A4I) in Appendix and K^ u (z) = K v {z) 
we finally obtain 



d 3 ' 2 



X 



V2V 5 / 4 A 3 r 



(33) 



Thus in the limit of narrow packets the transformed func- 
tion acquires a width of the order of A c 3> d. For all val- 
ues of d the width of the transformed function is larger 
than the width of the initial packet. 
For T z integrals in Eq. (|27p we have 



T z = V2(dV^) 3/2 
dT 



d 3 p 



e ipr/h e -p 2 d 2 /2h 2 cpz 



(2irh) 3 y / m 2 c i + c 2 p 2 



= -iK 



dz 



(34) 



Integrals T x and T y are obtained the same way. They 
have similar properties to T integral. 

It is seen that T of Eq. (33]) differs from D of Eq. (TJTD 
by a factor of (2y / nd) 3 ^ 2 . This is due to a different nor- 
malizations of the packet f(r) and the delta function. If 
one normalizes the packet according to J fs(r)d 3 r = 1 
the results of Eqs. ([35)1 and (|2T1) become the same. 

For the FW transformation, a function $(r) = 
/(r)(l,0,0,0) T is transformed to * = (S ,0,-S z ,-S x - 
iS y ) T , where So, S x , S y and S z are integrals defined be- 
low. There is [see Eq. (HU1) ] 



So 



d 3 p 



r/h fp(E p + mc 2 ) _ (d^) 3 ' 2 



(2-kK) 3 

[ ksin{kr)e~ k2d2/2 Jl + ^- dk, 
Jo V E k 



(35) 
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0.4- 



Gaussian 




FIG. 1: Elements in the first component of the transformed 
function *(r) (MOT) and $(r) (FWT) calculated for initial 
packet width d — X c . Dotted line — the initial Gaussian wave 
packet. 



where Ej- = + X^k 2 and f p is defined in Eq. ([T 
The integral S x is 



S x 



d 3 p 



^ipr/h 



fpCPx 



(36) 



and similarly for S y and S z . Integrals S x , S y and S z can 
be obtained as partial derivatives of an auxiliary integral 



d 3 p e ipr/hf p 

(2-Khf 



(37) 



with respect to x, y and z, respectively. 

In Figure[T]we plot various elements of the transformed 
wave functions: the integrals 2T (r), 2T z (z), S'o(r) 
and S z {z) calculated for d = X c . Solid lines repre- 
sent So(r) and S z (z) integrals, while dashed lines rep- 
resent 2T (r) and 2T z (z). The factor 2 in front of 2T (r) 
and 2T z (z) integrals is chosen to cancel out the factor 1/2 
appearing in Eq. (fT?]). The terms 2T z (z) and 2S z (z) are 
related to p z terms in Eqs. (JT7J) and (fH))) . 

The integral So(r) is very similar the initial Gaussian 
function. For the MOT the difference between the initial 
and the transformed functions is more pronounced. For 
larger d, the integrals So(r) and Tq(t) tend to the initial 
Gaussian function, so almost no widening occurs. This 
agrees with the the large-d approximation for Xb(r) in 
Eq. {2J|. 



IV. VARIANCES 

Finally, we estimate the spatial extent of functions sub- 
jected to the MO and FW transformations by calculating 



their normalized variances. For the MOT we define 

(ij-Vlf) 

VMO 



(38) 



and similarly for the FWT. Functions ^'(r) and $'(r) 
are normalized to unity. We assume the initial func- 
tion in the form of a Gaussian packet given in Eq. (|26p . 

Since *'(r) = *i(r) + %(r), see E q s - CH> and dHJ), 
it is convenient to introduce averages over the first and 



second parts of W(r). We define (r n ) a b 
with a, b= {1,2}. Then 



Vmo = 



(r 2 )ii + (r 2 ) 12 + (r 2 ) 21 + (r 2 ) 



22 



(r°)ii + (r") 12 + (r0) 21 + (r ) 2 2' 



(39) 



There is (r°) n = 1/2 and (r 2 ) n = 3d 2 /4. We consider 
now (r n )i2 terms. Introducing 



9 p n) = 



^ — ipr/H 



f(r)r n d 3 r 



(40) 



we have 



^•12 



d 3 P g P n) *f P ( 



t 9 

/ mc 



{2irh) 3 AE p 







icp z 



V 



icpr 
imc 2 + cp z 
cp + 



(41) 

The dagger denotes Hermitian conjugate of the vector 
and the star its complex conjugate. In the above equation 
the terms proportional to mc 2 cancel out. The terms 
including p z are odd functions of p z , so they vanish after 
the integration. For the same reasons there is (f* n ) 2 i = 0. 

To calculate (f") 22 we first introduce the Fourier trans- 
form of ^ 2 ( r ) 



2p 



e- lpr y' 2 (r)d 3 r = A. 



fjr, 



(42) 



where w p is the column in Eq. (1171) or the second column 
in Eq. gT]). Then 



J 4(2^)6 v 2pl v 2p 



A{2ttK) 
d 3 P |/p| 2 



4(2^)3 E 2 --^=2- (43) 

Since (r ) u = (r°) 22 = 1/2 and (r°) 12 = (r°) 21 = 0, the 
function ^'(r) is really normalized to unity. 
Now we calculate (r 2 )22 



f d 3 rd 3 pd 3 p' ir(r>- P ')/h 2yA & 
J 4(2^)6 V y? 2p'*2 P - 



(44) 



Integrating twice by parts over d 3 p and then integrating 
over d 3 r and d 3 p' one obtains 



{r )22 



= -h 2 



d 3 rd 3 pd 3 p' 
4(27rfi) 6 ' 

d 3 p ( fpW p 

A(2irh) 3 \ E p 



ir(p-p')/htf,A 



w 2p' v p*2p 



En 



(45) 



6 



After some manipulation we find 



(r 2 ) 



An 



where d = d/A c , and 



cxp(-i 2 d 2 ) t^dt 



(46) 



(47) 



For integer values of v and integer or half-integer val- 
ues of /j, the integrals A£ can be expressed in terms of 
modified Bessel and error functions, see Appendix. 
For the FW transformation we have similarly 



V, 



FW 



(<l|r 2 |<l) (r 2 ) FW 



where 

(r°) F w 



d 3 rd 3 pd 3 p' 



(48) 



(r°) F w ' 



(r )fw 



[ d 3 rd 3 pd 3 p> ir(p _ p0/n 2 - t ^ 
in which $„ = / e~ 4pr ' < I , ('r)<i 3 r. After some algebra we 



obtain 



' [A?_^-4A? /a ]. (51) 



In Figure [5] we plot calculated variances of the trans- 
formed functions divided by d 2 . For the initial Gaussian 
packet the variance is Vg = 3g? 2 /2 and it is indicated 
by the horizontal dotted line. The other two variances 
[MO-transformed function \P'(r) and FW-transformed 
function have very similar dependencies on the 

packet width d. Both variances have the zero-width 
limit Vd-j-o = 7d 2 /2 and they decrease with increasing d 
to the Gaussian limit Vg- 

Small-d and large-d limits of the variances can be ob- 
tained analytically using Eqs. (|46|) and ([51]) . Consider 
first the variance for (r 2 ) 22 from Eq. (|46[) . For small J the 
second term in Eq. (|46|) tends to zero (see Appendix) , and 
there is (r 2 ) 22 ~ lld 2 /4. The first term is (r 2 ) n = 3d 2 /4 
and the total variance is Vmo — (Hrf 2 /4 + 3d 2 /4) = 
7d 2 /2. Similarly, for the FW transformation at small d 
the second term in Eq. (|5"Tj) vanishes and the total vari- 
ance is again 7d 2 /2. 

The large- d limits of A£ in Eqs. (j46|) and (|5l"j) are cal- 
culated in Appendix. Applying these results to Eqs. (|46[) 
and (|5ip we find that the large-d limits of Vmo and Vfw 
are 3d 2 /2, i.e. they are equal to the variance of the ini- 
tial Gaussian function. Thus the widening of the trans- 
formed functions is pronounced for d < A c and it is small 
for d ^> A c , see Figure [2] It can be seen that, again, 
the FW transformation is more " compact " than the MO 
transformation because for a given packet width d there 
is always Vfw < Vmo- 
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FIG. 2: Normalized variances of the transformed functions 
for the MOT and FWT calculated for various initial packet 
width d. Solid line - FWT, dashed line - MOT. Dotted 
horizontal line — normalized variance of the initial Gaus- 
sian wave packet. The MOT and FWT variances has the 
small-d limit 7d 2 /2, while for large d they reach the Gaussian 
limit 3d 2 / 2 - 



V. CONCLUSIONS AND SUMMARY 



As mentioned in the Introduction, Rose Q demon- 
strated that a non-locality of the functional kernel gen- 
erated by the Foldy-Wouthuysen transformation is ex- 
tended in the coordinate space over the Compton wave- 
length A c . We generalize this result to the transformed 
functions and another transformation separating positive 
and negative energy states in the Dirac equation (Moss- 
Okninski transformation). In particular, we show that 
the delta function 5(r) is transformed into a nonlocal 
function smeared over the distance A c . Second-order 
variances are used to obtain similar results for narrow 
Gaussian wave packets indicating that the non-locality of 
the order of A c is a general property of the transformed 
functions. Our results strongly suggest that all "sepa- 
rating" transformations for the energies generate coordi- 
nate non-localities of this order. We emphasize that for 
relativistic-type equations, appearing either in narrow- 
gap semiconductors or in simulations, the non-localities 
can extend over hundreds of angstroms. 
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Appendix A 



We first quote formulas for the integrals appearing in 



the text. Let R = y/a 2 + b 2 . Then 



exp(itz)dt 
Vt 2 + a 2 



J a (bt)K a {ay/t 2 + i)tdt 



k s'm(kr) exp(—k 2 p 2 )dk 



2K Q {az), (Al) 
K X {R)/R, (A2) 



,-r7V 



(A3) 



x" : exp 



M 

-X- -7- 

4x 



dx 



4p3 

2(|)V_( M ),(A4) 



where 2?o(0> an d are the modified Bessel 

(MacDonald) functions. Next we collect formulas for A£ 
integrals. Let D = d 2 /2 and Erf(£) be the error function. 
Then 



A/2 = \eMD)K {D), 



A? = 



Al = 



7rexp((P)[l -Erf (J)] 
2 ' 
7rexp((P)(2(i 2 - l)[Erf(d) - 1] ^d 



(A5) 
(A6) 
(A7) 



For d — > there is dA®^ 2 ~ 0, because in this limit there 

is ^o(d) ~ ln(J) |2l|. Also, dA^ ~ and Jylg ~ 
because Erf(0) = 0. Using large- a; expansions for the 
Bessel and error functions we obtain 



lim 



A 



lim -^=A\ 



lim 



1 

r 
i 

2' 
1 

2' 



(A8) 
(A9) 
(A10) 
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